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Two kinds of configurations involving steps on surfaces are reviewed. The first one 
results from an initially planar vicinal surface, i.e. slightly deviating from a high-symmetry 
(001) or (111) orientation. In some cases, these surfaces separate into domains of different 
orientations by a mechanism which is very similar to phase separation in mixtures. The 
domain size initially increases with time and goes to a finite limit, whose value is related 
to elastic phenomena. The second kind of configurations results from making grooves in a 
high-symmetry surface. The surface smoothes out and takes an intermediate shape with 
facet-like hills and valleys, which are the source of a controversy which we try to clarify as 
much as possible. 



I. MOTIVATION. 

Surfaces are of great technological and fundamental interest since many processes, such as crystal growth, catalytic 
reactions, production of nanostructures, occur mainly at surfaces. 

The surface of an industrial catalyser is too complicated to allow for an observation of the microscopic processes. 
On the other hand, a high-symmetry, (001) or (111), surface of a cubic crystal is too simple. The only process which 
occurs is the diffusion of single atoms, which cannot easily be observed. 

The present review will be focussed on stepped surfaces of cubic crystals, i.e. surfaces made of large high-symmetry 
terraces separated by steps of thickness one monolayer (Fig. Q) . 

The fundamental interest of stepped surfaces is that they are both simple and not too simple. The ensemble of the 
steps can be arranged into many different structures, due to various factors which can play together or compete among 
themselves. Simple models can include attachment/detachment of atoms to/from the steps, diffusion or evaporation 
of atoms, and also step fluctuations. The steps can interact with each other, through short or long range interactions. 
The effect of a flux of atoms over the stepped surface can be considered, growth can occur through step flow, there 
can be step bunching, etc. The dynamics of steps gives informations on the microscopic parameters (energy, diffusion 
constants...) and are a test of our ability to solve problems of nonlinear dynamics. _ 

The experimental observation of steps is relatively easy, for instance by scanning tunnel microscopy (STM)EI. The 
present review will not be concerned with isolated steps, but only with experiments in which interactions between 
steps has an influence on their motion. Two examples are given in the next section. 
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II. EXPERIMENTS 



Solids are generally not in their equilibrum shape, and that is a good thing since otherwise one could not make, 
for instance, tools or artwork. However, small solids particles can reach their equilibrium shape at high enough 
temperature. 

A. Grooves 

A case of interest in the present review is that of a crystal surface in which periodic grooves have been made. Figure 
2 shows the example of periodically modulated Ni(llO) and Ni(100) surfaces which have been annealed at 800 C to 
obtain a steady state profile shape. This shape is sinusoidal in the first case but trapezoidal in the second. If annealed 
for long times, both profiles will decay and approach a flat surface. If the profile has a short enough wavelength and 
the temperature is high enough (wavelength A ~ 10/j.to , typical temperatures T between 900 — 1300-RT) the profile 
smoothens after reasonably short timesatl (Fig. ^|). 

The relevant questions in context with this experiment are: i) How does the amplitude 2h(t) decay? Exponentially 
or not?Ji) What is the transient shape of the profile? iii)How does the relaxation time r depend on the wavelength A? 
Mulling has given a simple theory (reviewed in chapter |HJ) valid above the roughening transition (we shall recall later 
what this means). Mullins' answers to the above questions, under the assumption of small slopes, are the following: 
i) The decay is exponential, ii) The profile is sinusoidal or becomes rapidly sinusoidal if it is initially not. iii) The 
relaxation time r is proportional to a power of A which depends on the type of dynamics. Mullins' theory will be 



reviewed in chapter III 



Early experimental observationscl evidenced that, for initially prepared sinusoidal profiles on Ni(100), (110) and 
(111) single crystals, the evolution towards equilibrium was found to be different depending on which surface the 
profile is created. For (110) surfaces the profile maintains the sinusoidal shape during the decay and the time law is 
strictly exponential. The relaxation time scales with the fourth power of the wavelength of the profile. In contrast, 
in the case of the high symmetry (111) and (001) surfaces, flat regions, reminiscent of facets, form at the maxima 
and minima. The decay of these profileswith time is slower than for the ones on Ni(110)En3. Non-sinusoidal profile 
shapes were also observed] for the (100) and (111) surfaces of gold shown in Figure 3. This case, however, is even 
more complicated than Ni, beca use vicinal Au(100) and (111) surfaces are unstable and facet. This problem will be 



discussed in parts [I B and XIII 



B. Phase separation on vicinal surfaces 

A vicinal or miscut surface is a surface obtained by cutting a crystal along a plane making a small angle a with a 
low index plane, (001) or (111) for instance. The ideal structure of a vicinal surface then consists of high symmetry 
terraces separated by equally spaced steps of atomic height. Indeed, such configurations have been observed, for 
example on vicinal copper surfaces at room temperatureau, or on vicinal Si(lll) surfaces at elevated temperature!!!! 
Although thermal fluctuations have been studied in detail on vicinal surfaces, and are responsible for a terrace width 
distributional, the average distance between steps is usually consistent with the ideal expected value. 

However, it has been known for a long timeO that upon thermal treatment, some stepped surfaces are unstable and 
separate into parts (or 'facets') of two or three different orientations. This separation is analogous to phase separation 
of unstable mixtures and the-,resultant morphology is a hill-and- valley structure. Such a situation occurs for instance 
on Au(lll) vicinal surfacesEj. Two examples are shown in figures |^ and [| namely Au(ll,9,9) and Au(5,5,4). Both 
surfaces are tilted by about 6° with respect to the (111) orientation and can be figured out as consisting of [1,-1,0] 
monoatomic steps which are initially equally spaced by about 24A. Both surfaces are unstable. After a thermal 
treament in ultra high vacuum, a faceted morphology is clearly seen in figures |] and g. In the simplest case, the 
unstable surface evolves towards an array of step bunches separated by high symmetry terraces. For instance, a 
(11,9,9) gold surface (Figure ^) transforms into (111) terraces of 35 Awidth, separated by small portions of another 
vicinal surface, namely the (7,5,5) orientation. In the case of the (5,5,4) surface (Figure^), the transformation leads 
to facets which are both portions of vicinals. All steps have still the same sign. The phase with the greatest slope 
corresponds to a (3,3,2) orientation, whereas the other one is misoriented from (111) by about 4°. The additionnal 
long range order displayed by the faceted morphologies will be discussed in section XII. The stability of the facets 
may be interpreted as caused by an interaction between steps which is not uniformly repulsive, as will be seen in 



subsection VB 



Another possible cause of faceting is surface meltingpJ'E-l For instance, Frenken et al.Ej have shown, using medium 
energy ion scattering, that vicinal Pb(lll) can separate into two surfaces of non-zero step density, a melted one and 
a dry one. _ 

Measurements^ of the faceting of vicinal Si(lll) surfaces show that this surface also rearranges into a hill-and- 
valley structure. At high temperature this surface consisted of an uniform density of steps and atJow temperatures 
two "phases" appear: the (7 x 7) reconstructed (111) facets and the unreconstructed step bunchesa. 

As in pcdijiary phase separation, the two coexisting phases can be deduced by the so-called double tangent con- 
structiontZHla if one knows the thermodynamic potential as a function of the relevant variable (the concentration in 
the case of a mixture, orientation in the case of a surface) . 

This thermodynamic potential, in the case of a stepped surface, can be in principle deduced, or at least understood, 
from the knowledge of the interactions between steps, which are the subject of section [v| 



III. VARIOUS TYPES OF KINETICS AND MULLINS' THEORY 



In order to model the experiments reported in Subsection [IB, we consider a crystal surface which is approximately 
a plane parallel to the xy direction, but for a small initial perturbation which slowly decays in time. We wish to study 
this decay. We can consider for instance an initially sinusoidal (unidirectional) profile of wavelength A, namely: 

z(x,t = 0) = fr(0)cos— — . (1) 
A 

The important quantity which drives the relaxation process is the chemical potential /i, i.e. the thcrmodyij 
potential per particle. If the surface is out of equilibrium, p is not constant on the surface and can be relates 
the surface free energy density by the Herring- Mullins formula. This formula can be written in a simple wayt3 if one 
introduces the free energy density tp per unit area of the projection onto any reference plane, e.g. the xy plane: 



p = p -a 2 



d ( dtp \ d f dip 
dx \dz x ) dy \dz y 



(2) 



where z is the height of the surface above the xy plane and z a — dz/dx a . This formula makes sense if tp is an 
analytical function of z x and z y . If it is so, the surface is called non-singular. The remainder of this section is devoted 
to the simple case of a surface which is non-singular in all its orientations. As will be seen in the next section this 
is not the case of stepped surfaces, but the forthcoming study will provide us the opportunity to introduce several 
useful concepts. 

For not too high temperatures, atoms can neither evaporate nor diffuse through the solid, so that surface relaxation 
takes place through diffusion of atoms on the surface. This case is called surface diffusion dynamics. One can then 
write a conservation equation which expresses the absence of evaporation: 

dz/dt= -a 2 V • j s . (3) 

where a is,the atomic distance. 

Mullinsu wrote the surface current density j s in (|J) as 

Js = -a- 2 f3D s V(i (4) 

where fi is the local chemical potential, f3 is l/(fcsT) and V = (d x ,d y ). The surface "mass diffusion constant" D s is 
defined as D s = a? p$D, where po is the adatom density at equilibrium and D is the adatom diffusion constant (we 
suppose that the advacancy contribution to diffusion is negligible compared to the adatom contribution). 

Mullinsu approximated ip by a quadratic function which will be assumed isotropic for the sake of simplicity. This 
procedure is usual in liquids and yields: 

c 2 a 

p{z x ,z y ) ~ — (z 2 x + z 2 y ) (5) 

where a is the surface stiffness0El and c is the interatomic distance along the z axis. 
Relations (0) and (§) yield 



/i = (i - a 2 c 2 a (z xx + z yy ) 



(6) 



Inserting (^|) into (Q) and (0), one finds 



dz/dt = -a 2 K V 2 (V 2 z) , (7) 
where _ftT = c 2 (5D s a > 0. This equation of motion is linear, so that an initially sinusoidal profile remains sinusoidal: 

z(x, f) = h(t) cos — - — . (8) 
A 

Using (0), the amplitude h(t) in (0) is seen to decay exponentially with the time t, 

h(t) = h(0)exp(-t/T) (9) 

with a relaxation time 

For an arbitrary, but periodic initial profile, the shape will anyway become sinusoidal after some time since it can 
be regarded as a sum of many sinusoids of different wavelengths (Fourier decomposition). At long times, only the 
component of longest wavelength persists. 

The above results hold for surface-diffusion kinetics. The comparison with experiment allows it to determine 



systematically the diffusion constant D on the (1 10) face of fee metals 

Another case is when matter is carried through the vapour. This process, important at high temperature, is called 
evaporation-condensation kinetics. For simplicity, the chemical potential of the vapour is always uniform and equal 
to the value /x of the bulk solid. The evaporation rate at a point of the surface where the chemical potential is fi is 
proportional to the difference /i — fj,Q : 

ioc^-^o (11) 
Then using expression @ for the chemical potential, one obtains 

Z jA {^Z XX -\- Zyy} (1^) 

where A is a temperature-dependent coefficient. Since ( |l2| ) is a linear equation, there are again sinusoidal solutions 
which satisfy eqs. (||) and (^. However, for evaporation-condensation kinetics, the above equations yield r cx A 2 in 
contrast with (|l0|). 

A third way to transport matter- and thus to allow for the relaxation of the profile, is bulk diffusion. When this 
process dominates, it can be shownta that, in the non-singular case, r oc A 3 . 



IV. STEPPED SURFACES 



The free energy of a stepped surface is the sum of the free energy of the high symmetry parts and the free energy 
of the steps. The simplest case is that of a regular array of straight, parallel steps at distance I (Fig. ||). If the step 
energy per unit length is 7 and if t is so large that the interaction between steps can be neglected, the projected free 
energy density resulting from high symmetry parts is independent of i while the free energy of the steps is "f/£. The 
total projected free energy density can be written as 

V (z x ,z y ) = V (0,0)+ 1 (z 2 x + z 2 y ) 1 / 2 (13) 

In contrast with ([5]), this expression is not analytic in z x and z y (Fig. ^b). A high symmetry surface is singular and 
Mullins' theory is not applicable. .— . 

More precisely, a high symmetry surface is singular below its roughening transitiom3 which occurs at a temperature 
Tr. Above Tr, 7 vanishes, ( |l3| ) does not apply, and one can argue that (|J) applies. For the (001) and (111) faces of 
most fee materials, Tr is very close or equal to the melting temperature. Low symmetry surfaces also have roughening 
transitions, which occur at much lower temperatures. 

In the case of stepped surface, the three types of kinetics (surface-diffusion, evapsijation-condensation, bulk diffusion) 
are still possible. However, surface-diffusion kinetics subdivide into two classesESc-a . i) In the first class, atoms are 
easily trapped by steps as soon as they meet one. The slow process is migration on the terraces. This is migration- 
limited surface-diffusion- kinetics, ii) In the second class, diffusing adatoms hardly feel steps. The slow process is 



attachment on steps (and detachment from steps). This is attachment-detachment-limited surfacc-diffusion-kinetics. 
We shall generally consider migration-limited surface-diffusion-kinetics. 

Attachment-detachment-limited surface-diffusion-kinetics is somewhat analogous to evaporation-condensation ki- 
netics, since the adatoms may be considered to have an uniform chemical potential fiQ as if there were a vapour 
bathing the surface. The exchange of matter between this mobile layer and the steps is assumed proportional to the 
difference between and the local chemical potential fi determined by surface energy changes resulting from surface 
displacements. \~\r~\ 

The equilibrium shape of crystals have facets parallel to singular orientatiorfiZrlla. As seen in the previous section, 
facet-like flat regions also appear in the transient shape of unstable grooves. A priori, both phenomena are expected 
to have the same origin. 

However, there are significant differences. At equilibrium, the chemical potential (i is the same everywhere. The 
size of a facet cannot decrease without locally increasing the density of adatoms (and therefore the chemical potential) 
near the facet, and the extra adatoms cannot be removed unless they go back to the facet, which thus comes back 
to its previous shape. On the other hand, in the case of the surface perturbations (not at equilibrium), the chemical 
potential deviates from its equilibrium value by a positive amount in the convex regions (top) and by a negative 
amount in the concave parts (bottom) and the removal of adatoms can be accomplished by transfering them to a 
bottom facet. Thus, if facets are for instance artificially made, they may be unstable with respect to exchange of 
atoms from top facets to bottom facets. Whether this occurs or not is a matter of dynamics: if the exchange of 
atoms between facets requires too much time, facets can persist until complete smooting. This dynamical problem is 
controversial, as wilLbe seen. 

A simple methodi^l is to replace the true free energy (|l3|)_( actually the improved expression ([lj) below, which takes 
step interactions into account) by a regularized form (Fig. |gbL-.This approach yields facets, and actually the transient 
groove shape is not very different from the equilibrium shapec3. However, the regularization procedure does not give 
much insight into what really happens near facets, and it is desirable to devise more microscopic methods. Before 
doing that, it is necessary to discuss interactions between steps. 

V. INTERACTIONS BETWEEN STEPS 

Formula ( |l3| ) consists of just the first two terms of an expansion. The next terms depend on interactions between 
steps. 

Interactions may be short ranged or long ranged. The definition which will be accepted in this review is that an 
interaction is called 'short ranged' if it decays faster than any power r~ a with the average distance r. Short range 
interactions generally depend on the microscopic features of each particular material and each particular orientation, 
while long range interactions have rather general features which make them especially interesting. 

A. Long range interactions 

Electrostatic interaction 

The most evident interaction between steps may be the electrostatic one: symmetry considerations suggest, tha t 
steps carry an electric dipole momentO. It is in fact a corollary of the electrostatic interaction between adatomsE§E2l. 
Actually, even a flat surface carries an electric moment density: it is the correction to this electric moment due to the 
step that we are considering here. 

The interaction between two dipole moments at distance r on the surface of a solid is proportional to 1/r 3 = 
(x 2 + y 2 )~ 3 / 2 , just as in a homogeneous medium. However, the proportionality factor is not easy to. evaluate. In the 
case of a metal, different evaluations have been given by Kohn & Lav£3 and by Lannoo & Allana. Note that two 
dipoles at the interface between two metals would have no long range interaction because of screening. 

If we come back to the case of steps, we can obtain the interaction energy between two straight, parallel steps at a 
distance £ by integrating the dipole-dipole interaction (£ 2 + y 2 )~ 3 / 2 over y. The resulting interaction is proportional 
to l/£ 2 , and of course to the length of the steps. The free energy density of a regular array of steps at distance £ 
is obtained by multiplying by the step density which is itself proportional to l/£ = (z 2 + z 2 ) 1 / 2 . The free energy 
density per unit projected area ip (introduced in formula^) is then, for \z x \, \z y \ <C 1 

<p(z x , Zy) = <p(0, 0) + 9l (z 2 x + z 2 ) 1 / 2 + g 3 (z 2 x + zlf' 2 (14) 



where g% = 7 is the step free energy per unit length and the constant g% = <?f can be positive or negative, depending 
of the orientation of the dipoles. 



Elastic interaction 

Another type of long range interaction arises from the elastic stress which is unavoidably exerted by a step (like any 
other defect) on the surface. This stress-pipduces a strain which is of course maximum near the step, and decreases 
as 1/r 2 with the distance r to this stepE2iliZlii3 . As a result, an elastic interaction energy between two parallel steps at 
distance i arises, which is-jrapportional to l/( 2 , just as the electrostatic interaction! For steps of identical sign, this 
interaction is repulsivcHEZlllj. When it is taken into account, the resulting surface free energy density of a vicinal 
surface has still the form ([l4]), but 53 now contains an elastic contribution gf as , which is always positive in contrast 
with gf ec . |_. 

A noteworthy exceptior£2l to the l/£ 2 rule is the elastic interaction between steps in the case of an adsorbed 
epitaxial layer which has a different lattice parameter from that of the substrate. In this case, the (elastic) dipole- 
dipole interaction has to be integrated over the whole misfitting terrace, not only over the step edge: the interaction 
becomes thus logarithmic! 



Entropic interaction 

A third type of interaction between steps arises from the fact that steps cannot cross each other. This can be 
regarded as an infinite contact repulsion. If the free energy density is calculated as a function of the average distance 
I between steps, this contactr*epulsion turns out to generate a term proportional to £~ 3 , just as in the case of elastic 
and electrostatic interactionsE3~E3. This term may be understood as resulting from the reduction of the entropy of a 
step by its neighbours. Indeed, the entropy is essentially the logarithm of the number of degrees of freedom and a step 
looses degrees of freedom every time it meets a colleague, since it has no right to cross it. This effective interaction 
may therefore be called entropic. It is obviously repulsive since each step tries to have as much space as possible. 

When the contact repulsion is taken into account, the projected free energy density ip has S ^rtih^|fi^p*|-(H), but 
a new term gf™*, which is positive as g^ las , has to be added to obtain g$. The calculations^SSEaa'tsTof <?J"* 
uses the analogy between an array of mutually-avoiding steps and the trajectories of non-interaction, one-dimensional 
fcrmions. The resultEa is 

ent K 2 {k B T) 2 

9s =—^- (15) 

where 7 is the line stiffness. This result is twice as large as that obtained by Jayaprakash et a@ because of a factor 
1/2 in the second term of equations (3) and (13) of these authors. 

All three interactions of g% = gf ec + gf as + <?f nt are equivalent to a force proportional to l~ 2 between pairs of steps, 
or to a term proportional to in the free energy density. However, there are important differences: the entropic 
interaction is a repulsion between nearest neighbour steps only while elastic and electrostatic interactions relate all 
pairs of neighbours. 



Experiments 

The case of the solid-supcrfluid interface of helium is particularly accessible to experiments because equilibrium is 
much more rapidly established than in other materials. Certain experimepta^j suggest an interaction which decreases 
slower than r~ 2 . A similar suggestion has been made-by Arenhold et afl for Pb. A theocetical attempt to explain 
that have been made (in the case of He) by Uwahael On the other hand, Balibar et alcil have pointed out that 
experimental results must be very cautiously interpreted. 



B. Oscillating and short range interactions 

As said before, short range interactions are generally particular to each material and it is difficult to say anything 
general. Oscillating interactions have similar effects since they have local minima at short distances. One can mention 
the following mechanisms. 



Oscillating electronic interactions 



For certain metal surfaces there is a partially unfilled, electronic surface band. Then a phenomenon analogous to 
Friedel oscillations generates an oscillating interactioncj between point defects. This interaction is proportional to 
r~ 2 cos(2fcp<r), where &f is the Fermi wavevector and r the distance between defects. Integrating along a row generates 
an interaction between parallel line defects (e.g. steps) which is proportional to r^ 1 cos(2kpr). Since it oscillates with 
distance, the absolute minimum corresponds to a finite value of the distance Iq between steps. The effect is very 
important since surface orientations with steps larger than £q can become unstable as observed experimentally in 
certain cases (see subsection II B). Oscillations are washed out when r becomes longer than the electron mean free 
path, and therefore the interactions between parallel steps on a clean surface are likely to decay as 1/r 2 at long 
distances in all cases. An oscillating interaction appears also without partially unfilled electronic surface band, but it 
decreases faster with distancec 2 ! 



Interactions resulting from reconstruction 

Surface reconstruction may have a similar effect, i.e. make certain orientations more favourable than the interme- 
diate one, in other words, create an oscillating interaction between stepsJ13 

A typical example of reconstruction is provided by the (001) surface of silicon or any semiconductor. Semiconductors 
usually have the diamond structure, in which atoms have a very low coordination (4 neighbours). Creating a (001) 
face in a rigid lattice would reduce the coordination of surface atoms to two neighbours and that would cost a lot 
of energy. Therefore the surface atoms form 'dimers', thus reducing the symmetry. This symmetry reduction with 
respect to that of the rigid lattice is called surface reconstruction. Now, step positions which would cut dimers are 
unfavourable, and this disavantages certain distances between steps. I— . 

The orientations which are favoured by reconstruction are sometimes called 'magica. 



VI. MEAN FIELD THEORY FOR A WELL-CUT SINGULAR SURFACE 



After having recalled the various general properties of surfaces and steps, we come back to the experiments on 



groove smoothing on Ni(001) and Ni(lll) described in subsection II A . The case of Au(001) and Au(lll), whose 

vicinals are unstable, is excluded here and will be addressed in section XIII 

Steps are complicated objects and in this section we shall yield to the temptation to consider them as straight 
lines (Fig. [7] a). This is only possible if the average surface is initially perfectly cut along a symmetry orientation 
and if there are no dislocations. Even then, steps are not straight, but fluctuating lines (Fig. [?] b). Neglecting-these 
fluctuations may be_called a "mean field approximation" . This approach was chosen by Rettori & Villained and 
Ozdemir & Zangwilo, who assumed an interaction energy proportional to L/£ 2 between nearest neighbour steps of 



length L at distance I. These features are those of the.£ntropic interaction addressedin subsection V A. 

We will recall here the results of Rettori and Villained and Ozdemir and ZangwilEl Below Tr, the analytic form 
(||) of the projected free energy <p is not correct. Instead, one should use (O), where 33 is given by ( |l5| ) in the absence 
of step interactions except contact interaction. In a continuum approximation, the instantaneous shape of the surface 
can be described by a height function z(x,t) (if the perturbation is unidirectional) and ( |l4| ) takes the form 

<ys = 9l\Zx\ +93\zx\ 3 (16) 

Using (§) and (§J), one obtains the excess chemical potential asEJEj 

8fi = -6a 2 g 3 \z x \z xx (17) 

where z xx = d 2 z/dx 2 . 

If the relaxation occurs by surface diffusion, we can still write (0) and (Q), and together with ( p7| ) the equation of 
motion is obtained 

Yt = -G^W*- ( 18 ) 



where G3 cx g^D 'g/ksT and D s = o?pqD as before. Eq. (|l8|) is non-linear and until now no analytic solution has been 
found for it. I— . 

Since <5/i should vary smoothly in space, ( |l7| ) is consistentcZI with a singularity near the top (Sx ~ 0) of the type 



z = h{t) - C{t)\Sx\ 3 / 2 



(19) 



where h(t) is the amplitude of the profile at time t. 

The singularity (|f9|) is the same as predicted at equilibrium near a facet. However, according to Rettori & Villain 
and to Ozdenir & Zangwill no facets are expected to form during flattening. Indeed, the decay of the top (and bottom) 
terraces (should they form facets or not) is due to the junction of the top (or bottom, respectively) ledges. When two 
top (or bottom) ledges approach from each other, why should lower ledges stay behind instead of keeping contact with 
them? If top ledges keep contact with lower ledges, no facet is expected in the decaying profile. Therefore, according 
to (|l9|), the profile is predicted to be sharper than a parabola at its top (and at its bottom as well). 

Ozdemir & Zangwill have solved (|l8|) , or more precisely a discrete version of ([l8]) , and found under certain assump- 
tions (shape preserving solution) that the amplitude h(t) decays with time t as 

h(t) = h(0)—^— (20) 

1 + t/Tx 

where t\ ~ A 5 . Relation ( pp| ) is often quotedEfi but it is not always mentioned that t\ has to depend on h(0). 
Indeed, if ( |20| ) holds for a particular value hg of h(0), with a particular value r" of t\, it follows, for any time t\ 

(o < h < t), 



l + ii/r° l + t/T x yu l + t/r° x 1 + (t - h)/(h + r°) 1 + h(h)(t - h)/(h T°) 

If the time origin is translated to t\, this coincides with (p0|), but t\ is proportional to l//i(0): 

rx = r° x h o /h(0) (21) 

The weakness of the above theory is that any interaction between steps of opposite sign (at the top and the bottom 
of the profile) is ignored. This would be acceptable if the ledges were straight (Fig. | a). In reality, they have 
fluctuations (Fig. R b). We shall try to take these fluctuations into account in section Ixl. 



VII. FACET-PREDICTING THEORIES OR THE AMBIGUITY OF THE CHEMICAL POTENTIAL 

A paradox of the mean field theory outlined in the previous section is the energy release which results from the 
pairwise annihilation of ledges at the top and bottom of the profile when they are assumed to be straight. Spohnc3 and 
Hager & SpohnE^I have proposed a completely different theory of groove smoothing which does not show this paradox, 
at least in the case of surface diffusion kinetics. In the case of evaporation-condensation kinetics, both approaches 
predicts the same singularity, a sharpening of the profile at maxima and minima for long times : z ~ |<5x| 4 ' 3 (this 
prediction is only valid in the limit of a slowly varying surface profile, i.e., a minimum sample size is needed). 

The remainder of this section is devoted to surface diffusion kinetics. In this case, Spohn's results regarding the 
shape of the profile are in strong contrast with those described in section [Vl|. He finds, indeed, that facets form in 
the transient regime. ■— ■ 

As pointed out by Tangpa, the reason of the disagreement lies in different expressions of the local chemical potential 
/iona surface which is not at equilibrium. Since the current is a function of \7fj, as seen above, the consequences are 
crucial. 

The chemical potential \i is the sum of its value in the solid far from the surface, and an excess S/i resulting 
from the surface tension. 

If a system has N atoms and a free energy F, its chemical potential is /i = F/N. 

Now, let us consider the topmost terrace of the profile of Figure [?] a. If its width is £ and if the sample length is 
L, then the number of atoms in the terrace is N = iL/a 2 while, if interactions between steps are neglected, the free 
energy excess resulting from the surface energy is SF — 2jL where 7 is the step free energy per unit length. If the 
whole terrace is removed, the corresponding variation of free energy is Nfio + SF hence S/j, = SF/N or 

5fx = 2^a 2 /e (22) 

Experiments are performed on macroscopic samples, where this expression is much larger than the term coming 
from step interactions. 

If ( |22| ) is applied to grooves, the same formula holds for the bottom terrace with an opposite sign, 5 fx = —2"fa 2 /£ 
if the bottom terrace has the same width. The average chemical potential gradient is therefore V/i = Aja 2 /(£X) and 
the current from the top to the bottom is of order 4D70 2 / (£X). If £ is small, the bottom terraces are rapidly filled and 



the top terraces are rapidly peeled, and after some time, only the intermediate terraces are left, and facets appear at 
the top and at the bottom. This is Spohn's description. 



In section VI, on the other hand, an underlying idea is the following: if the edges of the topmost terrace move so 
that the width changes by small amount di, the number of atoms of that terrace changes by dN = 2Ld£/a 2 , but 
the total length of the edges does not change, it is always 2L, and therefore the chemical potential is, disregarding 
interactions between steps, Sfi = 0/dN or 

5fj, = (23) 
Instead of a straight step, it is of interest to consider also a closed, e.g. circular terrace of size R. It contains 



N = kR /a , where a is the atomic area, and its free energy excess is 5F = 2njR. The formula analogous to (22) is 

Sfi = 6F/N = 2-ya?/R (24) 
while, arguing along the lines of section VI, one obtainsEll 



2 

= d(5F)/dN = (25) 
R 



Formulae (|24| ) and (25) are still different, but just by a factor 2. The results which can be deduced from both 
formulae are just quantitatively different, not qualitatively. In any case smaller terraces have the higher chemical 
potential and therefore decay first. 

Formulae (|24|) and ( p5| ) are in contradiction. Which one is correct at equilibrium? It depends. 

For a single terrace on a high symmetry surfac e, (B5|) is correct for a circular terrace of radius R, while 5/i = for 
a stripe-shaped terrace, as can be obtained from (|25|)by replacing the ledge radius R by oo. 

However, if a finite crystal has an approximately circular facet of radius R at equilibrium, S/i can be obtained by 
removing a whole atomic layer from a facet. The step interaction energy is not much modified in that operation, and 
therefore (^4|). If the crystal is obliged to have the shape of a bar, S/i is given at equilibrium by (p^). Of course, it is 
also given by (|l7|), and both expressions should be equal at equilibrium. 



Out of equilibrium, there is no reason a priori to write ( |22j ) or (|23|), (|24j) or (25). A careful analysis of the dynamics 
is necessary to solve the controversy. 

Relation (p5[), which relates the chemical potential to the curvature, is a two-dimensional form of the Gibbs-Thomson 
formula. 

An objection to Spohn's theory is that it does not describe how pairs of top (or bottom) steps merge. The peeling 
of extreme layers are accompanied by fluctuations which are not better taken into account by him than by Rettori & 
Villain and Ozdemir & Zangwill. In the next section we shall try to find a solution of the controversy in simulations. 



VIII. SIMULATIONS 

The relaxation of profiles have been simulated in (1+1)D and (2+l)D by many authors&@'0El. The solid-on-solid 
(SOS) model have been used by most of them. Its Hamiltonian is 

U = J l^-^WI ( 26 ) 

xy,x'y' 

where h xy denotes the height of the surface at site xy, the sum runs over nearest-neighbours pairs of sites and J is 
the bonding energy. 

Both diffusion and evaporation dynamics were considered in the simulations, and in general the standard Metropolis 
algorith was used. 

For evaporation-condensation, a conserved dynamics is sometimes used, in the form of a detailed balance of particle 
exchange between the surface and the surrounding gas: a surface atom can move from its initial position to some 
other, arbitrarily distant site (conserved number of surface atoms). In some works, non-conserved dynamics are 
used. For (2D4JJ and T > T^, the simulation results agree well with the predictions of Mullins: the profile relaxes 
as a sinusoicOtilJ. For T < Tr it was^bservecEl that the profile fluctuates, in time, between a sinusoidal and 
a "trapezoidal" form. Selke and Biekero observed that these fluctuations are accompanied by two distinct time 
scales in the flattening procedure, reflecting the step wandering and the shrinking of islands. They assumed that the 
meandering stage leads to the dominant time scale and their results for the relaxation time r are compatible with 



t oc A 3 , but they speculate that for larger systems sizes, r oc A 4 . Actually, recent large scale simulations by TangE2l 
confirmed the exponent 4, but with a correction in the power law. 

For diffusion, the number of particles must be conserved and papticlfi exchange occur between nearest-neighbors 
sites. For temperatures_above Tr the results of the simulationsEJCSO agree well with the predictions p£ .Muffins 
(sinusoidal profile decayed, the valueja = 4 was obtained for the scaling law r oc A™ of the relaxation timeE3'L3). For 
temperatures T < Tr, some authorsEj found that the scaling behaviour breaks down after some time and then cusps 
and plateaus develop at the top and bottom of the profile, othcrsElo found that the shape of the relaxation profile 
fluctuates in time between having a flat top to having a rounded top (sinusoidal). We remark that these authors have 
considered typically small initial profile heights h — 4, 5. Different scaling laws for the typicakrelaxation time r were 
obtained: n = 4 (obtained by Searson et alE3), n = 3.5 — 4 (obtained by Erlebacher et alO), n = 5 (obtained by 
Ramana Murty et alo). 

We mention that most authors performed their simulations using standard Monte Carlo algorithms, and from that 
they obtained time dependent laws for the quantities of interest. Presently it becomes more cleaioEl that the kinetics 
of the relaxation must be preferable studied by methods such as the Kinetic Monte Carlo (KMC), which includes a 
prescription to take into account the real time instead of to take the unit aLtima as the interval between two updatings, 
for instance. Only very recently KMC simulations have been performedEjE^H on the relaxation of grooves. 

Another important point is that, since in the evaporation-condensation dynamics atoms are exchanged between 
the surface and the vapour and in diffusion dynamics it is necessary to consider the motion of the atoms from a site 
to another (in general a nearest-neighbour), surface diffusion takes much more computing time than evaporation- 
condensation. Due to that, in many simulations (for diffusion dynamics) some short-circuit artefact have been 
introduced, especially, small samples have been considered (small heights, or small widths). 

In what concern the choice of small samples, one of the problems is that to obtain the correct scaling laws (e.g. the 
relaxation time r vs the wavelength A), the amplitude h and the average distance between steps i = a\/(Ah) must 
be large with respect to the atomic distance a. 

What is a short sample? 

The problem related to the width L of the system in the direction of the grooves is different. In available analytical or 
numerical theories, elastic or electrostatic interactions between the steps which form the p rofile are not considei^-ajid. 
then the only interaction which remains is the contact repulsion. As seen in section V A, in-analytic theoriesOEJO'LJ 



this contact repulsion is taken into account by introducing an effective free energy equal mm i 



to 



fent-L^ (27) 

for each pair of steps at average distance £. 

Here, 7 is the line stiffness related to the energy W per kink by the relation 

Since steps are not straight at finite temperatures, it is useful to define the average distance £ between kinks along 
a step, given by 

tK .( 1+ sSm). (29) 



Formula (|27]) is correct only for long distances £. Moreover (27) makes sense-only for large L. An exact calculationEa 
-for instance considering a single fluctuating line between two straight linea£3 at distance 21— shows that ( p7f ) holds 
only if the relation 

L > ii 2 /a 2 (30) 

is satisfied. This yields a first condition on £ 

i < L(4h/Xf. (31) 

Whether this condition is satisfied or not depends on the temperature T, which should anyway be smaller than 
Tr. An order of magnitude of Tr is given by writing that the free energy of an isolated step per atom, which is 
W — fcsTlnfl + 2 exp(— (3W)] at low temperature, vanishes. Thus 



exp — — ~ 2 



(32) 



As a matter of fact, if one wishes to check the analytic predictions made for singular surfaces, it is safer to choose 
T much lower than Tr. Indeed the aspect of a surface on short distances is very similar just below Tr and just 
above Tr. In particular, closed terraces are present between steps, and steps exhibit "overhangs" which are generally 
ignored in available theoretical treatments. The choice T < Tr/2 seems reasonable. This implies, according to ( |29| ) 
and (p^), a second condition on £ 

£ > 3a (33) 



It turns out that conditions— ( and ( p3| ) are not always simultaneously satisfied in simulations, specially when 
diffusion kinetics is consideredSH The simulation by Ramana Murty and CooperEJ is a special case, since they used 
sufficient long samples, and actually Ramana Murty and Cooper do find the scaling law r oc A 5 predicted by Ozdemir 
and Zangwill. Note, however, the very low value of h in all simulations. The choice of all lengths (A, h, L) results 
from a compromise between too large values which would saturate computers, and too short sizes which would have 
no relation with reality. ._. 

A conclusion which can be drawn from simulations^ is that facets do appear at least in certain cases, but it is 
not quite sure that this phenomenon is not related to the small sample size. For a given sample length, facets seem 
to be observed in simulations at low temperature, when_the steps have but few contact points. On the other hand, 
at higher temperatures, the profile is roughly sinusoidaS. In the following section, an explanation of the results on 
short samples is presented. 



IX. THE TRANSIENT ATTRACTION 
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The interactions between steps we considered so far are thermodynamic. They correspond to terms of the free energy, 
and have an effect on the equilibrium state. When writing kinetic equations, these thermodynamic interactions appear 
in the detailed balance relation 

pi eM-PE A ) = Vb eM~f3E B ) (34) 

where denotes transition probability per unit time from a state A to state B and Ea is the energy of state A. 

In a non-equilibrium situation, we will prove that a clustering of steps of identical sign can occur in the case of 
surface diffusion kinetics. This can be interpreted as arising from an attraction between such steps, even though there 
is no thermodynamic interaction, i.e. Ea = Eb has the same value for all configurations of interest. 

attraction can only be done for a 1-dimensional surface (Fig. a) which 



Unfortunately, a quantitative theory of t 
is a good representation of a short sample. 

The one-dimensional model is fully characterized by the probabilities per unit time and unit length, a + (£) and 
a~(£), for two neighbouring steps at distance £, to exchange an atom in one direction or in the other one. The case 
of interest is when both steps have identical sign. Consider the configuration B obtained from a configuration A by 
transfering one atom downward from a step to the neighbouring one at distance I. 

The transition probability per unit time from A to B is, by definition, a + (£). The transition probability per unit 
time from B to A is aT (£ + 2a) because the terrace width in B is £ + 2a. In the absence of interaction between steps, 
the detailed balance relation writes 

a-(£ + 2a) = a+{£) (35) 

If the distance I between both steps is large, most of the emitted atoms go back to the step where they started 
from. The transfer probabilities <X—l£) and a + (£) are therefore expected to decrease with increasing t. 
More precisely, it can be showrH, that: 



^ = ^lTUVa (36) 



and 



where ao is a constant and the "Schwoebel length"E-3i23 £ s (equal to in most of available simulations) depends on 
the detachment probability of atoms from steps. The emission rate aoa per lattice site is obtained if one notices that, 
at equilibrium, adatom emission by steps is compensated by absorption. The absorption rate is proportional to the 



product of the adatom density po by the adatom jumping rate 4D/a 2 , where D is the adatom diffusion constant. 
Therefore 

a = . (38) 

a 

There is no thermodynamic interaction between steps at equilibrium since ( |36| ) and ( |37| ) satisfy the detailed balance 
relation (|35|). However, a~(£) > a + (£), which means that atom exchange between two steps of identical sign tends to 
shorten their distance. It is hardly avoidable to interpret this as a kind of attraction between neighbouring steps of 
identical sign. 

Of course, all steps should be taken into account, and the effect disappears for a regular array of equidistant steps 
(an equilibrium configuration) but it is present for the periodic profile of Fig. || a, as demonstrated by Monte-Carlo 
simulations^! The result is a blunting of the profile; steps of identical sign attract themselves, so that the maximal 
slope increases, and consequently the top and the bottom flatten. 

The following exercise helps to understand the nature of the interaction. Consider a one-dimensional profile made of 
alternating pairs of up and down steps (Fig || b) . Let the distance between two consecutive steps be initially uniform 
and equal to £q. It will first be assumed that steps of opposite sign cannot annihilate -an unphysical model which 
has the advantage to have a non-trivial equilibrium state in which the average distance between two consecutive steps 
is £q. In this model the average distance < £(t) > betwen two consecutive steps of identical sign will first decrease 
because of the attractive interaction, until collisions take place. Then, as an effect of the contact repulsion, < £(t) > 
will increase and come back to its equilibrium value £q. Thus, the attraction has just a transient effect. For this 
reason, we will call it transient attraction. 

In a one-dimensional profile with alternating groups of 2h/c up steps and 2h/c down steps, the transient attraction 
does produce an initial clustering of identical steps, which produces facets, which are observed in Monte-Carlo sim- 
ulationsH. In these simulations, steps of opposite sign are allowed to annihilate, and it turns out that facets persist 
until the profile completely flattens. This is a surprise, but the effect seems to be observed only on one-dimensional 
samples or in two-dimensional ones which are too short in the groove direction. Therefore, the transient attraction 
will be ignored in the following sections. 



X. STEP FLUCTUATIONS AT THE TOP AND BOTTOM 

Analytic theories of Rettori and VillainEl and Ozdemir and ZangwillEl are mean field ones, where fluctuations are 
neglected. Away from the top and bottom, step fluctuations are expected to be reasonably taken into account by 
the entropic interaction. The challenge is now to introduce step fluctuations of the top and bottom ledges into these 
theories. The following calculatioriEil is an attempt to solve this problem. 

Because of these fluctuations, finite terraces, i.e. closed step loops, form at the top of the profile, which have a 
typical length L\ (Fig. ^b). The width of the loop should be comparable with the distance 2£\ between both topmost 
infinite steps of each period. Each finite terrace terminates by a tip whose radius is of order £-y, so that the excess 
chemical potential before the tip is 

S^=^f (39) 

The velocity v of the tip is proportional to the current density of adatoms flowing from the tip which is equal to 
(3poD\7fj,. Since the tip is at distance £\ of the terrace edges, the order of magnitude of Vp, should be expression (|3£ 
divided by £\ , so that 



PpoDja 4 



(40) 



£\ 

The time t\ necessary to peel the upper layer is obviously 

ri = Li/v (41) 



so that (|o|) yields 



The nucleation time of a finite terrace after the previous one has disappeared is expectedEl to have the same order 
of magnitude r as the lifetime of a finite terrace. 

Topmost steps can only fluctuate by emission of atoms. These atoms can only go to the lowest steps. This suggests 
that step fluctuations arise from the exchange of atoms between topmost and lowest terraces. This exchange is of 
course biased, and preferably from the top to the bottom. But it is also partly stochastic. Its mean square fluctuation 
will be calculated without taking the bias into account. The exchange rate in each direction is expected to be 
given by a formula analogous to (|36|) or p7[ ) where the distance I is replaced by the half- wavelength A/2 while the 
proportionality factor is given by (|38|). Thus, the average number dv / (dtdL) of atoms exchanged from the top to the 
bottom per unit length in the unit time is (neglecting bias) 

dv _ 8p D ^ 



dtdL A 

The average < v > of the number v of particles emitted by a length b of a topmost step in the time t\ is the 
product of this quantity by br\. The mean square fluctuation < Sv 2 > of this number is, if the various absorption 
and emission processes are independent: 

< Sv 2 >=< v >= -^r&Ti (44) 

This fluctuation of the number of atoms produces a fluctuation 81 of the distance between the topmost steps on a 
length b. Expression ( f44| ) can therefore be identified with the square of the area bS£, so that 

..o a 4 < v > dv Tia 4 
<S£ > -^ 2 — = dML— (45) 



The length b should be chosen as the most efficient one, which maximises (45). Therefore, b should be chosen as 
short as possible. However, expression ( fl5"| ) cannot be larger than the thermal fluctuation, so that < 8£ 2 >< kgTb/j. 
The appropriate choice of b is therefore 

t < S£ 2 > 

*>= J h T (46) 
k B T 

The topmost steps meet when < 5i 2 >= I 2 . This relation, together with ([t5]), ([42]), @) and (|||), yields an 
expression for L\\ 

-y 2 A£ 2 

L > " 5P* < 47 > 

The average current density j from a top terrace is obtained by calculating the total current from each tip, which 
is proportional to the product of ( |4C| ) by I, and dividing by the distance L\ 

(48) 

This can be identified with the current density of atoms from the top to the bottom terrace, which can also be 
roughly evaluated by dividing the chemical potential difference 25 /i (given by formulae [jj] and |l9|) by A/2 (in order 
to obtain the chemical potential gradient) and then multiplying by the same kinetic coefficient p$(3D as in fi3). 

. = 27C 2 p pDg 3 a 2 (4Q) 
2A 



Identification of fl48|) and (|49|) yields 



^ = (50) 



27C 2 g 3 £ 



l 



Eliminating Li between (pOj) and (E^), one obtains 



276^337 



If <?3 is replaced by its value (p"5|), this relation reduces to £ 3 » C~ 2 . This is just what could be deduced from ( |l9| ) 
by replacing h — z by 1 and &c by £1. In other words, the results of the mean field theory of section VI would not be 
considerably modified if fluctuations were taken into account. 

An approximate expression of C(t) as a function of the modulation amplitude h(t) can be obtained if one assumes 
that (19) is approximately correct for z = 0. Then, since A 3> i\, 



C(t) « h(t) (A/4)- 3/2 

The peeling time of the top layers is given by inserting (pf!) and (pi) into (f42|). One obtains 



(52) 



Tl 
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8fc|T 2 f3p D-fa 4 8k 2 B T 2 (3 Po Dja 



and, using 



16(fc B T) 2 



4/3 



16(fc B T) 2 



27^37 



ft(t) 



-8/3 



4/3 



0-8/3 



(53) 



The quantity which can easily be measured is the time after which the amplitude h has been reduced by a factor 2 
for instance. This decay time is 



T de 



(54) 



The power A 5 is the same as obtained by Ozdemir and Zangwill (see formula |2l]) but expression ( pi"| ) was proportional 
to The amplitudes h used in Monte-Carlo simulations are too small to decide which formula is correct, and so 
far as we know, there is no experiment which can prove or disprove (p4j). 



XI. MISCUT SURFACES 



Metal surfaces cannot easily be cut along a high-symmetry direction. There is generally a small miscut angle a. If x 
is the direction common to the high symmetry plane and to the miscut average surface, the groove direction is defined 
by its angle (Vr/2 + ip) with x. A sinusoidal surface modulation causes the steps to assume a specific shape. For if,' = 
they are sinusoidal in an projected on-top view while for ip = it/2 they are straight, with their separation modulated 
sinusoidally. Staying with a low amplitude/wavelength ratio (surface slopes small compared to a) and using a surface 
free energy of the kind in eq.(|T^) solutions for the rate of profile decay via migration-limited, surface diffusion kinetics 
were derived by Bonzel and MullinsEa for two particular cases: ip = and ip = n/2. In the first case the rate is 
proportional to (g\ + 3<73tan 2 a)/ sin a, which is approximately equal to <7i/sina because g% tan 2 a <C gi- Thus the 
step self-energy is the important energetic quantity, driving the profile decay. This is understandable because the 
wavy steps want to reduce energy by becoming straight which is equivalent to a reduction in profile amplitude. For 
ip = tt/2, on the other hand, the decay rate is proportional to 653 sina/ cos 4 a, i.e. it is driven by the step interaction 
energy only. Since the steps are basically straight from the beginning, the total energy can only be reduced by a 
redistribution of all steps, so that their separation becomes about equal. This redistribution of steps causes the 
amplitude to decrease. The ratio of rates for these two types of modulation is approximately 51/(6(73 sin 2 a) for 
a < 10° assuming equal surface diffusion coefficients normal and parallel to steps. Since g\ is believed to be larger 
than <?3, this ratio can be a large number, depending on the choice of a. An experimental verification of the effect was 
given for 1.5° and 5° miscut Au(lll) surfaces that had been modulated in the ip = and ip = n/2 modes. Figure || 
shows the experimental amplitude decay at 1023 K of two profiles in orthogonal orientations on the 5° miscut Au(lll) 
surface. The substantial difference in decay rates illustrates the predominant influence of step and step interaction 
energy, respectively. 

Strictly speaking, since there are unstable orient ations near the (111) facet of Au, the theory by Bonzel and Mullins 
is not valid. This case will be addressed in section XIII. 



Unfortunately, this kind of experiment always yields products of surface diffusion coefficient times surface energy 
term, and unravelling these two factors seems difficult. A combination of kinetic and equilibrium experiments offers 
a possibility to separate them but the total effort is substantial. I— . 

Profile decay is exponential with time on vicinal surfaces for the small slope approximationta, and the relaxation 
time r is proportional to A 4 . This is the same result @ as for a non-singular surface. 



A numerical investigation beyond the small slope approximation has also been dondrj for the ijj — n/2 case. The 
calculated decay is linear with time. It also does no longer obey a 1/sina dependence but rather saturates with 
decreasing miscut angle (Fig. |l0|). For tp — 0, the corresponding profile shape is sinusoidal for small slopes and 
more trapezoidal when larger slopes are permitted. This blunting can be understood intuitively as follows o Near 
an isolated step, 5[i is proportional to the step curvature as seen from (p5[). Therefore, the steps try to decrease 
Sfi by increasing their radius of curvature, where it is most easy, namely at the top and at the bottom, while in 
the intermediate region they are squeezed together, so that they cannot change their shape much. The increase of 
the radius of curvature results in a blunting of the profile. For small amplitudes, the squeezing effect is absent and 
the profile remains sinusoidal. Another interpretation of the blunting is the following. The miscut has the effect to 
suppress the singularity of the free energy as a function of z x ([l6|), so that the Bonzel-Preussca calculation applies, 
which does predict a blunting. Simulations showed that increasing the step interaction energy relative to the step 
energy causes the steep portion between maximum and minimum of the profile to become less steep, as one expects. 

Keeping the step energy parameters constant, the influence of a on profile shape was also studied using the same 
algorithm. The smaller the miscut a, the flatter the tops and bottoms of the profile. The behavior of the interesting 
tp = tt/2 case for larger slopes, where steps of opposite sign would be generated by the modulation, has not been 
dealt with so far. The question whether such profiles would facet or not remains unresolved. This case is, from an 
experimental point of view, an enormous challenge because the positioning of the profile parallel to the intrinsic steps 
would have to be vepjj accurate. 

Bonzel & MullinsE§ have also treated the decay of a periodic profile of small amplitude (h/X <C a) when there is 
a crossover from migration-limited to attachment-detachment surface-diffusion kinetics. The crossover occurs in the 
vicinity of a wavelength A*. The above relation r cx A 4 , characteristic of migration- limited kinetics, holds only for 
A 3> A*, while r cx A 2 , characteristic of attachment-detachment-limited kinetics, holds only for A « A*. This is a fairly 
remarkable result which contrasts with the.-jcase of competition between surface diffusion kinetics and evaporation- 
deposition kinetics on non-singular surfaces!! Then, r is proportional to A 4 for short wavelengths and to A 2 for long 
wavelengths. 



XII. LONG RANGE INTERACTIONS 



To our knowledge, all computer .simulations have been made without taking long range interactions into account, 
while theoretical treatmentsc3c3'HE3 use an entropic interaction in l/£ 2 between neighbouring steps, which corre- 
sponds to a contact repulsion. _ . 

This is a questionable choice, because MetoisEi has experimentally shown that the elastic or electrostatic interaction 
is larger than the entropic interaction even at fairly high temperature. 

In the present section, we briefly discuss the expected effect of long range, repulsive interactions... and find quite 
disappointing results. 

As seen in section V A , the total elastic and electrostatic interaction per unit length between two steps at distance 
£ is C/£ 2 for steps of identical sign, and C 1 jl 2 for steps of opposite sign, where C and C are constants. 

We shall first adress the case C — C > 0, neglecting step fluctuations. Then, the sign of the steps is irrelevant, 
and therefore a state with equally distant steps (Fig. [ll]) is an equilibrium state, since it is an equilibrium state when 
all steps have the same sign. For appropriate step signs, one can give to this equilibrium state any arbitrary height 
h and wavelength A. The existence of an equilibrium state of arbitrary height and wavelength in the absence of step 
fluctuations suggests a very slow decay when step fluctuations are taken into account. The profile has a zigzag shape 
without facets. This picture is very far from experiments. 

Can we do better by assuming C > C > 0? It is easily seen from a perturbative treatment that there is still 
an equilibrium state of arbitrary height and wavelength in the absence of step fluctuations, at least if C is not too 
different from C . The profile has again a zigzag shape without facets. The tops and bottoms are even sharper. We 
do not get closer to experiments although the fluctuation-induced relaxation can be somewhat faster. 

In the case C > C > 0, the perturbative treatment shows that there is also an equilibrium state of arbitrary height, 
but now the profile has facets. However, an extremely slow relaxation is expected. 

These results are quite disappointing. On one hand, steps should interact. On the other hand, taking interactions 
into account yields results which disagree with experiments. This paradox disappears if the surface is miscut. 



XIII. UNSTABLE ORIENTATIONS 



In section 0, we mentioned that vicinal (001) and (111) faces of Au are thermodynamically unstable. This property 
may be related to an interaction between steps which is not uniformly repulsive in contrast with the assump tions 



made, for instance in section VI . The origin of this interaction may be one of those addressed in subsection |V B . 
However, although it is well known that Au(lll) possesses a 22 x V3 reconstruction, the role played by the surface 
reconstruction is not so clear as in the,case of vicinal Au(100) surfaces or vicinal Si(lll) surfaces. "Magic vicinals" 
have been invoked, both theoretically^ and experimentally, in order to explain the remarkable stability of particular 
vicinals. Such surfaces are characterised .by a strong relationship betweepa_.the terrace width and the reconstruction 
mesh. It has been observed for exampleOlJ on vicinals Au(100) and alscOll3 on vicinal Pt(100). The role of surface 
reconstruction is clearly evidenced in the pioneering work of Williams et al.EJ on Si(lll), which is reported in subsection 



IIB. 



As pointed out initially by MarchenkcO and more explicitly by E. Wi lliam s and her colleaguesEJ, elastic effects are 
important in the "phase separation" phenomenon described in subsection II B . Without taking elasticity into account, 
an unstable surface, e.g. Au(ll 9 9) would separate into domains whose size would increase indefinitely with time. 

Elasticity limits the size of the domains, because each-surface element of each "phase" exerts a stress on the solid, 
the result of which is an elastic interaction proportionalo to l/|r — r*| 3 between two surface elements d 2 r and d 2 r' at 
points r and r*. The total energy of a stripe of size £ is obtained by summing over the two coordinates of the points 
atoms r and r". This operation multiplies by the sample width L y and introduces 3 integrals which transform the 
function l/\r— r*| 3 into a divergent function of £ proportional to hi£. Similarly, the energy of a periodic array of stripes 
of two different orientations is proportional to L y ln£ per period, so that the elastic energy density is proportional 
to hx£/£. The coefficient turns out to be positive if the average orientation is unstableOEl. Since the energy density 
resulting from domain boundaries is proportional to l/£, too large domains are thermodynamically unstable with 
respect to domain wall formation. To summarize, an unstable surface can reach an equilibrium state by forming a 
periodic structure of alternate "phases" . The value of the period is given by a balance in energy between the elastic 
energy and the energy cost of bou ndar ies. Indeed, a long range order is present on the stepped surfaces of Au(lll) 
described previously in subsection [IB (figures || and ||). The morphology of the faceted Au(ll,9,9) surface (figure 
^) is clearly periodic. The value of the period is about 65 A, as it was determined both using diffraction technique 
(LEED) and statistical analysis of STM images. In the case of Au(5,5,4), the quasi-period is about 1800 A. The factor 
of more than 20 between the values of both surfaces can be qualitatively understood as follows. 

The surfaces of interest are vicinals of (111) defined by Miller indices (mmp), with m/p close to 1. They all 
correspond to [1,-1,0] steps, but these steps have a different structure for m/p < 1 and for m/p > 1 (compare fig.^a 
and |b). Their energy can therefore be quite different. Moreover, since the elastic energy is proportional to hi£, the 
equilibrium domain width £ eq is an exponential function of the step energy, andrJthis can explain the very high ratio 
between the widths in figures fig||a and [|b. Details are given by Rousset et al.LJ It should be pointed out that, so 
far, no evolution of the faceted morphologies has been found by increasing, either the temperature of annealing (since 
400°C until 700°C), or the time of annealing (since a few minutes until 24 hours). This suggests that the observed 
morphologies are close to the equilibrium state. Perhaps not quite the equilibrium state. Indeed, the evolution toward 
the thermodynamically stable state can be expected to be very slow. This slow evolution results from the attractive 
interaction between steps which has been seen to be presumably responsible for the instability. This interaction 
favours the formation of step bunches observed in Fig. |4| The motion of these bunches is presumably more difficult 
than that of isolated steps, just as big particles diffuse less easily in any medium than slow particles. Moreover, if 
the bunches are close to the optimal size £ eq , the evolution toward this size requires breaking the bunches, which 
is certainly very difficult. More precisely, the evolution toward equilibrium should be rather slow until the period 
reaches a value of order £ e q/2, and much slower yet afterwards. The above discussion holds for mmp surfaces with 
m/p < 1, whereas the case m/p > 1 is more complicated. i_. 

The attraction between steps is also expected to produce a very slow decay of grooves on Au(lll) and Au(001)E3. 
These grooves are completely different from those on Ni(lll) and Ni(001), whose vicinals are stable. In particular, 
there is no objection to the formation of facets since unstable orientations result frorn-attractive interactions between 
steps of identical sign, while the instability of facets, advocated by certain theoristsc3c3, is related to the repulsive 
interaction. Moreover, these facets are expected to decay very slowly since the upper and lower steps are tightly 
bound to their colleagues and can but hardly, detach from them and combine with the opposite step of opposite sign. 

Groove healing has also been observedlljO on Au(lll) vicinal surfaces (a = 1.5°, 5°), where the decay of pertur- 
bations of orientations ip = and ip — ir/2 were observed (ip is defined in section XI). These profiles could not be 
treated under the small slope approximation. Initially the etched profiles have nearly square wave shapes, but after 
some time of annealing, sharp edges are eliminated. The decay of the "0 = profiles is much faster than the one for 
ip = tt/2 profiles. The former shows a well-rounded shape after some time of annealing, while the ip — n/2 profile 



does not become sinusoidal even after a very long time. 

For ip — ir/2 a (111) facet is present, in adition to an extra facet which is flat bounded by edge-like features 
indicative of missing orientations. The profile in this case becomes very asymmetrical. A possible explanation is that, 
when the profile is etched on a vicinal surface, the terrace widths are not the same. When the terrace widths are 
smaller than a given minimum, if the average orientation of this part of the profile is close to a magic orientation, it 
is energetically more favorable for the surface to present a flat region than to maintain the structure of steps/small 
terraces. 



XIV. CONCLUSION 



We have reviewed the state of the art of groove smoothing. 

When the average surface is non-singular and far from any singular orientation, everything is clear and there is 
nothing to add to Mullins' theory. 

When the average surface is singular and the neighbouring orientations are unstable, the situation is pretty clear 
too. Facets should appear at the top and at the bottom, they have sharp edgesj23'EJ and the decay is expected to be 
generally very slow. An experimental check of the theory is not yet available, and can obviously be made only if the 
decay is not too slow. 

When the average surface is singular and the neighbouring orientations are stable, e.g. Ni(OOl), the situation is 
more complicated, except for a miscut surface . 

Let the controversial state of the art be summarized in the case of a singular surface of an ideal, perfectly cut 
crystal, when the neighbouring orientations are stable. The most common theoretical approach is a generalization of 
Mullins' theory, using the non-equilibrium chemical potential fi. However, the choice of the appropriate formula for /u, 
is controversial. A correct theoretical treatment should correctly take into account the peeling of the successive layers, 
and therefore the fluctuations of the top and bottom ledges. The attempt presented in section [x| tends to confirm 
the results obtained without taking fluctuations into account, i.e. there is no facet on a well-cut profile. However, 
the treatment of fluctuations in section [x| is far from being exact. The most receni experiments suggest that facet 
do appear in well-cut surfaces as well. This is in agreement with certain theoriescSO, although these theories do not 
take fluctuations into account in a satisfactory way either. 

A point which requires further investigation is the 'transient' attraction investigated in section IX. It is clearly 
effective for an (unphysical) one-dimensional surface. It appears between two steps of identical sign because an atom 
diffusing from the lower to the upper step has a shorter path and therefore a lower probability to come back. It does 
produce a blunting of a one-dimensional profile. Whether it produces a similar effect on real surfaces is still unclear. 

Most of existing theoretical treatments of the groove problem ignore long range (elastic or electric) interactions 
between steps. The attempt which was done here does not yield satisfactory results. However, these interactions 
are extremely important in the other problem reviewed here, namely phase separation of vicinal surfaces. A careful 
comparison of both types of experiments in the same materials would probably lead to a better understanding of 
interactions between steps. 
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FIG. 1. Stepped areas found on a Au(lll) surface misoriented by 50'. Note that steps are straight on the right and lower 
part of the figure, whereas kinks are present on the left and upper part. All steps are monoatomic, 2.35Ahigh. The size of the 
STM image is 1020Ax 1020A. Light white lines, running perpendicular to the steps, are the sign of the 22 x %/3 reconstruction. 



FIG. 2. Interference micrographs of periodic profiles on Ni(110) and Ni(100) single crystal surfaces, both annealed at 1173 K 
for several hours. The profiles exhibit a wavelength of 12 micrometers and amplitudes of 0.21 and 0.27 micrometer, respectively. 
Note the completely different shapes in the two cases. From Yamashita et al. 



FIG. 3. Interference micrograph of a periodic profile on a Au(lll) surface annealed at 1123 K. The wavelength and 
amplitude are 7.0 and 0.125 micrometer, respectively. 



FIG. 4. Decomposition of an unstable surface: Au(ll,9,9). (a) Hard sphere model of the ideal Au(ll,9,9). Note that the 
steps consist of {100} microfacets, displayed in grey, (b) STM image of the Au(ll,9,9) faceted surface. The size of the image 
is 356Ax 314A. All steps are monoatomic (2.35A high) and aligned along the zone axis < 0, 1, —1 >. 



FIG. 5. Decomposition of an unstable surface: Au(5,5,4). a) Hard sphere model of the ideal Au(5,5,4). Note that the steps 
consist of {111} microfacets, displayed in grey, (b) STM image of the Au(5,5,4) faceted surface. The size of the image is 1.07 
fi m x 1.07 fi m. (c) zoom in of the frontier between both phases, and located at the white spot on fig. 5b. On the left part, 
terraces are 13A wide, whereas on the right part they are about 35A wide. 



FIG. 6. a) A vicinal surface, b) Free energy as a function of the slope 



FIG. 7. a) Decay of grooves considered as arrays of straight steps, b) Effect of fluctuations 



FIG. 8. a) The one-dimensional surface of a two-dimensional crystal, b) The same as in the previous figure when the height 
is just two layers 



FIG. 9. Log-plot of profile amplitude versus annealing time for two profiles of 4.3 micrometer etched on a 5° miscut Au(lll) 
crystal in the ijj = iv/2 and if) = orientations. The ratio of decay rates at 1023 K is 25. From Surnev et al. 



FIG. 10. Numerically calculated decay rate versus amplitude of periodic profiles with modulation if) = for three different 
miscut angles a. Parameters: step free energy = 0.4, step interaction energy = 0.16, wavelength = 4 . From Bonzel & Mullins 



FIG. 11. A stable configuration (in the absence of fluctuations) in the case of steps interacting through a long range interaction 
independent of the sign of the steps. 
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